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Abstract 



We propose a density functional for anisotropic fluids of hard body particles. It interpolates between 
the well-established geometrically based Rosenfeld functional for hard spheres and the Onsager 
functional for elongated rods. We test the new approach by calculating the location of the the 
nematic-isotropic transition in systems of hard spherocylinders and hard ellipsoids. The results are 
compared with existing simulation data. Our functional predicts the location of the transition much 
more accurately than the Onsager functional, and almost as good as the theory by Parsons and Lee. 
We argue that it might be suited to study inhomogeneous systems. 



1 Introduction 

The density functional approach is one of the most powerful and widely applicable ap- 
proaches to nonuniform fluids pL Its idea is to express the free energy as a functional of 
locally varying one-particle densities F{p{r)}, with an ideal gas contribution and an excess 
free energy which accounts for the interparticle interactions. This allows one to calculate 
the structure and properties of fluids with various kinds of inhomogeneities. Density func- 
tionals for simple fluids have reached a considerable degree of sophistication and are able to 
describe such complex phenomena as freezing, wetting, and surface melting 

One of the most established density functionals for hard sphere fluids is the Rosenfeld 
functional Q. It reproduces by construction the Percus-Yevick solution for the correlation 
function, which is known to be very good ||. As a consequence, it provides highly accurate 
predictions for the structure of inhomogeneous hard sphere fluids [|]], It has been used suc- 
cessfully to study mixtures of hard spheres J5|] and polydisperse fluids |J. After a simple 
anisotropic mapping procedure it can be applied to fluids of fully aligned ellipsoids 0]. Ex- 
ploiting the Gauss-Bonnet theorem, it has been generalized for isotropic molecular fluids [||]. 
This generalized functional has been used as a starting point to calculate direct correlation 
functions of isotropic multicomponent fluids §. A Rosenfeld type approach has been de- 
veloped for mixtures of rods and needles, assuming that the needles are too thin to interact 



with each other directly [10], and for systems containing only one single rod immersed in a 



fluid of spheres [11]. However, to our knowledge, the Rosenfeld functional has not yet been 



extended to general anisotropic fluids. 



The simplest density functional for anisotropic particles is the Onsager functional [12], 
which truncates the virial expansion after the leading coefficient. Onsager showed that 
this functional produces a nematic-isotropic phase transition []l^] in fluids of sufficiently 
elongated particles. Its predictions are in good agreement with experimental results on 



systems of tobacco mosaic viruses [13]. One can even show that it describes the transition 



exactly in homogeneous systems of infinitely elongated rods |1J]. Thus any functional for 



nematic liquid crystals should reduce to the Onsager functional in this limit. 

In the present paper, we propose a density functional for hard anisotropic particles which 
interpolates between the Rosenfeld functional and the Onsager functional. It reduces to the 
Rosenfeld functional in the case of particles with spherical symmetry, and to the Onsager 
functional for homogeneous fluids of infinitely elongated particles at density close to zero. 
As a first test of the functional, we have calculated the nematic-isotropic transition for 
hard spherocylinders and hard ellipsoids, and obtained reasonable results. We believe that 
our functional might provide a useful new approach to the study of inhomogeneous liquid 
crystals, e.g. the study of interfacial phenomena such as surface anchoring. 



2 Background 

We consider a fluid of hard anisotropic particles with positions r and orientations w. The 
grand canonical free energy is the minimum of the free energy functional 

(3n[p) = (3F d [p\ + (3T ex [p] + J drdw pit, w) [n - V(r, w)] (1) 

with respect to the one-particle density p(r, w), where (3 = is the Boltzmann factor with 
the temperature T, p the chemical potential and V(r, w) summarises all external potentials. 
The first term describes the ideal gas contribution 

pjM = J dr dw p ( r) ro ) ^ ln[p(r, w)X^\ - l) (2) 

with the de Broglie wavelength At- The second term T ex is the excess free energy, the 
central quantity in density functional theories [jj]]. 

In the Onsager theory, the excess free energy of two hard particles is given by 

pfff = - - J drdw dr' dw' p(r, w) p(r', w') M(r - r ; , w, w'). (3) 

The Mayer-function M takes the value (—1) if two particles at (r, w) and (r', w') overlap and 
vanishes otherwise. In the homogeneous case, the one-particle density p does not depend 
on r and can be written as p(r, w) = p f(zu). The expression then reduces to 



N 2 



J dw dw' f(w) f(w') v(zu,m'), (4) 



where u(w,w') is the covolume of two particles with orientations w and w' , and iV the 
total number of particles. Eqn. (||j) corresponds to a virial expansion up to second order. As 
mentioned in the introduction, the Onsager theory is exact in the limit of infinitely elongated 
particles |fM[| . Compared to real thermotropic liquid crystals, it tends to overestimate the 
shape anisotropy required to observe a nematic phase at a given finite density. 

A natural extension of the Onsager model would be to include higher order terms in the 
virial expansion. On principle, this is feasible, and an extension of the Onsager functional 
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up to third order was actually carried out for ellipsoids [15]. However , the calculations 
are very cumbersome due to the increasing complexity of the integrals. As long as one 
is interested in homogeneous systems, one way to overcome the problem is the decoupling 
approximation. It consists of a resummation of the virial expansion, where the first virial 
coefficients are calculated exactly, and the remaining ones are approximated by a mapping 
onto a reference system. Given a virial expansion of the form: 



oo l 

£ P -B l+1 [f (w)] , 

1=1 L 



(5) 



the decoupling approximation reads ||16| - |19| |: 



If 



£ P -B l+l [f H] + 



i=i 
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(6) 



where Bi+i [f (w)] are the virial coefficients of the anisotropic fluids, calculated exactly up 
to the r + 1 order and ref refers to the reference system. Two possible reference systems 
have been proposed: the hard sphere fluid and the isotropic fluid of the particles under 
investigation. If r = 1 and the reference system is the hard sphere fluid with the Carnahan- 
Starling equation of state pofl , we recover the Parsons-Lee functional [16, 17] 



PL 



1 (4-3*;) 
8 (1 -r/) 2 



dr dzu dr' dm' p(r, w) p(r', w) M(r — r', w, w'), 



(7) 



where rj is the packing fraction. The Parsons-Lee functional predicts accurately the location 



of the nematic-isotropic phase transition for a wide range of shape anisotropies |21, 22 1. It 
can also be applied to inhomogeneous systems. However, it does not describe the microscopic 



structure very well: it yields very crude correlation functions [23, 24 1 and produces rather 



unrealistic density profiles in inhomogeneous systems as a consequence |25|[ . 

The success of the Parsons-Lee theory in predicting the isotropic-nematic coexistence den- 



sities motivated Poniewierski and Holyst |26| and Somoza and Tarazona |27j to combine the 



decoupling approximation with a "weighted density approximation" (WDA) scheme [28, 29, 
|30f] - Both are constructed such that one automatically recovers the Onsager functional in the 



low density limit. In the hard sphere limit, Somoza's and Tarazona's [27] and related [31] 
versions reduce to the well-known Tarazona functional |^9|| . The latter incorporates the 
Carnahan-Starling equation of state and the Percus-Yevick direct correlation function for 
homogeneous hard sphere fluids, and has been used with great success to study inhomo- 
geneous hard sphere fluids and solids. Applied to systems of spherocylinders, the exten- 
sions |2^, [3l| of the Somoza functional generate a nematic phase and several smectic phases. 
Graf and Lowen have put forward a simplified "modified weighted density approximation" 
(MWDA), and used it to reproduce complex phase diagrams of spherocylinder fluids 

Besides the Tarazona functional, another equally successful hard sphere functional has estab- 
lished itself in recent years: the fundamental measure theory or Rosenfeld functional 0, f|, g[. 
It has the advantage of being based on somewhat more fundamental considerations: it does 
not require the explicit input of the equation of state and the direct correlation functions, 
instead they pop out automatically. One obtains by construction the Percus-Yevick solution. 
Moreover, the functional is formulated a priori for general convex particles of arbitrary ge- 
ometry. It thus seems to call for a generalization to anisotropic fluids. Unfortunately, there 
is one problem: the original Rosenfeld functional does not reproduce the Onsager functional 
in the low density limit. 
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Here, we propose a simple modification of the Rosenfeld functional which solves that prob- 
lem for anisotropic particles, and reduces to the original Rosenfeld functional for isotropic 
particles. Like Somoza's and Tarazona's functional, our functional interpolates between an 
established hard sphere functional and the Onsager functional. Thus we hope that it will 
be equally useful. 

Before introducing our approach, we sketch briefly the basic equations of the original Rosen- 
feld functional ||. It can be formulated for general multicomponent fluids of convex hard 
particles. In our case, different components i may be identified with different particle ori- 
entations. The functional then reads: 

= jdr( $x(r) + <D 2 (r) + <& 3 (r) ), (8) 

*i(r) = -no(r)ln(l-n 3 (r)) (9) 
* / \ ni(r)n 2 (r) -ni(r) • n 2 (r) 

*»w = — ^) ( 10 ) 

$ , rl _ H(r)-n 2 (r)(n 2 (r)-n 2 (r)) 
3[) 8vr(l-n 3 (r))2 {L±) 

The $s depend on the weighted densities 

n«(r) = Y,J dr'pi(r) ^ a) (r - r ), (12) 

i 

with the number density of the ith component Pi{r) and the weight functions 

^ 0) (r) = a(r)5(R t (r)-r)K t (r)/4iT wf ] (r) = Q(r) 6(R t (r) - r) 

wf\r) = Ci(r) 5(Ri(r) - r) fli(f)/47r ^ 3) (r) = 9(^(r)-r) 

w 4 (1) (r) = Ci(r) 5(Ri(r) - r) n<(r) ##)/4tt wf(r) = Q(r) 5(Ri(r) - r) fi^r). 

(13) 

Here f denotes the unit vector in the direction of r, Ri(r) the radius from the centre of 
a particle of type i to the surface in the direction r, Hi(r) and -fQ(r) the local mean and 
Gaussian curvature on that specific surface point |34, [35|], n«(r) the outward unit normal 
(see Figure 1), and 5 and 6 are the usual delta and step function. The factor Cj(r) ensures 
that integrals over the weight functions wf 1 (v) and w£ (r) (a ^ 3) are surface integrals 
over the surface of the particle i. For spherical particles, it can be omitted. In general, it is 
given by 



C«(r) = ^det gi (r)/detg iiSph (f), (14) 

where gj(r) is the local metric tensor of the particle i, and S pj 1 (r) that of a reference sphere 
of radius Ri(r). (In polar coordinates, ^det gj S pi 1 (r) = Rf(r)s'm(9)). Finally, the factor A 

in Eqn. ( pi] ) ensures that the dimensional crossover for hard spheres is reproduced correctly 
and that the hard sphere system exhibits a solid- fluid transition in three dimensions 0|. 



WW 2 ) 3 ... , n 2 (r) 



A= l-3£(r)' ^ ^ 



n 2 {r) 



(15) 



Alternatively, Eqn. (11) can be replaced by a more modern version due to Tarazona [33], 
which depends on an additional tensorial weighted density and does not contain the ad hoc 
factor A. 
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Figure 1: Illustration of R(r) and n(r). 

3 Construction of the functional 

As mentioned before, the Rosenfeld functional describes hard sphere fluids very successfully. 
On principle, it could also be applied to molecular fluids. However, it has the serious 
drawback that it does not reduce to the Onsager functional in the low density limit. Even 
worse, a closer inspection of Eqns. (js|)-(|l~5|) reveals that it has no contribution at all which 
would favor parallel alignment of particles in a homogeneous fluid. Hence it cannot produce 
stable homogeneous anisotropic fluids. 

The reason for this failure can be understood by looking at the relation between the Mayer- 
function and the weight functions in more detail. For hard spheres, the Mayer-function for 
a pair of particles at position rj and rj can be decomposed exactly as Q 

i\/r (0)k> (3) , (3)^ (0) . (1)„ (2) . (2)„ (1) (1) _ (2) (2) _ (1) / 1e x 

— Mij = w\ ' ®w y j ' +w\ ' <g>u>j +w\ ; <g>w y j +w\ ' ®w)j -w\ y (g)wj ' — w> ; <g>wj ; , (16) 
where (8) denotes the convolution product: 



W\ ® Wj' 



[ drw (a \r-ri)w^\r-rj). (17) 



This decomposition together with Eqns. (|8|)-(|il^) ensures that the functional reproduces the 
correct virial expansion at least up to second order. Unfortunately, a decomposition of the 
type above is no longer valid for anisotropic particles. Therefore, the Onsager limit is not 
recovered. Chamoux and Perera (9| have proposed ways to cure the virial expansion on 
the level of the direct correlation function, for the case of isotropic molecular fluids A 
systematic way of dealing with the problem on the level of the functional itself would be to 



add a correction term on the right hand side of Eqn. fll6|) , to deconvolute it (if possible), and 



then rederive a Rosenfeld type functional on the basis of this new decomposition. This can 



be done for rod-sphere interactions An extension to the general case is currently under 



way |3q]. Unfortunately, it turns out that infinitely many additional terms are required in 



the decomposition in order to make Eqn. (16) exact, and the numerical treatment becomes 



difficult unless one resorts to approximations. 

Here we propose a simpler, numerically perhaps more tractable Ansatz - a straightforward 
modification of Eqns. (||), ( |i0| ) and (pTT|) . Instead of keeping n\ and ni as separate quantities, 
we suggest to replace the products rti(r)ri2(r) and ni(r)-ri2(r) by new joint quantities Bi2(r) 
and ni2(r) such that the functional reproduces correctly the Onsager limit. This is achieved 
as follows: we define the functions 

g(u7i,Wj,Ti - Tj) = gij := -My - ® u>j 3) - wf' <8> ivf 1 (18) 

h(muWj,Ti - rj) = hij := ® wf + wf ] <g> wf - wf ® wf - wf ® wf (19) 
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and 



H eS (wi,Wj,ri - rj) = gij/hij. (20) 



Then we replace ni(r)ri2(r) and ni(r) • n2(r) by 



If rr 

nu(r) ■= - dmid-cuj dridrj p(ri,uji) p(rj,Wj) H e (m^ zjj, - r^) (21) 
x (wf'fc - r) Wj(Tj - r) + wf'fa - r) wy(rj - rfj 

If rr 

"-12(f) := - Jdwidwj dxidij p{ri,wi) p(rj,Wj) H {za^w^Vi - rj) (22) 
x (wP (rj - r) • wj 2) (r,- - r) + wf } (r { - r) • (r,- - r)J . 



rr 

The effective function H is constructed such that the functional reduces to the Onsager 
limit (g) in the low density limit. Thus the mean curvature H(r) in the weight functions 
it/ 1 ) and w^ 1 ) (Eqn. (p^)) is effectively replaced by a function H ■ H e ^ which depends on 
the orientations of a pair of particles, and their distance vector. In the hard sphere limit, 
H is constant, H = 1, and we recover the original Rosenfeld functional. Note that n\i 
and h±2 contain information on pairs of interacting particles. We thus give up the idea of 
formulating a functional which depends only on single, orientation independent weighted 



densities [37]. An approach in the same spirit has been introduced by Schmidt for mixtures 



of spheres and needles ||10| . 



The functional can be simplified by making the approximation that H e ^ only depends on 
the orientations W{ and Wj. In that case we solely require that the second virial coefficient 
is exact in the homogeneous fluid for every distribution of orientations. In the homogeneous 
fluid the contributions to the free energy density having vectorial character vanish. Carrying 
out the integrals on the spatial variables, we obtain 

PFr _ Wl . , pRSJdw dw' f(w) f(w>) H eff ( ro , w>) 1 p 2 S 3 
N ~ [ V)+ 1- V + 24vr(l- ?? )2' [Z6) 

where R, S are the mean radius and the surface of the body, respectively, and 77 is the 
packing fraction. The mean radius is defined as the integral of the mean curvature over the 
surface of the particle, R = J s dA H/4tt. Eqn. ( p3|) must reduce to the Onsager functional 
([|) in the low density limit p —> 0. With rj = vop this leads to the equation 

H^^)= ^ W ^ /) - V0 , (24) 
v ' RS K ; 

which replaces Eqn. (|20|) in this approximation. 

In the case of homogeneous fluids, both (|^) or ( |24[ ) give the same bulk free energy as a 
function the orientation distribution f{w): 

^ - U1 n \ V i P fadwif{w)f(m')v(m,m') 1 p 2 S 3 

IT' V) l^ + 2 T^n + 247(1^)2 W 



Note that for isotropic fluids, f(w) is constant, and the free energy (25) is identical to that 
obtained with the original Rosenfeld functional. 
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Figure 2: Illustration of a spherocylinder 

4 Application to spherocylinders and ellipsoids 

We have tested our approach by calculating the location of the nematic-isotropic phase tran- 
sition for hard spherocylinders and uniaxial prolate ellipsoids. The results were compared 
with the corresponding phase diagrams obtained from the Onsager theory, the Parsons-Lee 



theory, and simulation data [£2], [3^, |39 |, 



A spherocylinder consists of a cylinder of length L and diameter D capped by a hemisphere 
of diameter D at both ends (see Figure 2). 

The excluded volume of two spherocylinders which have the angle 6 with respect to each 
other is 1 12] 

v{9) = 2DL 2 sin(0) + 2vr D 2 L + 4ttD 3 /3. (26) 

In the case of ellipsoids, the exact calculation of the excluded volume is quite involved. 
One possibility is to follow Ref. [fly], another to exploit the Perram and Wertheim routine 
for the ellipsoids contact function |3l]]. Here, we have adopted a scheme outlined in Ref. 
[22], because it is sensibly faster. We consider a pair of equal uniaxial prolate ellipsoids 



with semiaxes a, a and b; one ellipsoid is mapped onto a sphere of radius a, and the other 
one is mapped onto a particular biaxial ellipsoid which has always one semiaxis equal to 
a. The remaining two semiaxes depend on the relative orientation of the two ellipsoids and 
can easily be evaluated. The excluded volume between the sphere and the biaxial ellipsoids 
can be calculated using an expression by Kihara for the orientationally averaged excluded 
volume between two convex hard particles H3]. The covolume of the original pair of uniaxial 
prolate ellipsoids is then recovered by inverting the mapping. In the past, many theoretical 



studies niq, 17, 43, 14, ktq] have considered systems of "hard Gaussian overlap" particles, i. 



e., ellipsoids where the contact distance is approximated by an expression due to Berne and 



Pechukas |46]. However, it has been noted some years ago ||T^] and shown recently [47] that 



the comparison of these theoretical results with true ellipsoids behaviour is not appropriate. 

Both in the cases of hard spherocylinders and hard ellipsoids, the properties of the fluids 
are fully determined by the density and the shape anisotropy parameters, x = L/D and 
k = b/a, respectively. 

In the homogeneous fluid, the orientation dependent part of the excess free energy has the 
same form in the Onsager model, the Parsons-Lee model, and our model: 



(3^[f{m)] _ A 
N 2 



J dm dm' f(m) f(m') v(m,m'), (27) 



where A = p in the Onsager model, A = f (i-^)a m the Parsons-Lee model, and A = 2(1-77) 
in our model. Once a model free energy functional is chosen, the next step is to determine 
the thermodynamically stable phase for a range of densities. We thus need to calculate the 
functions f(m) that extremize the free energy under the constraint jdm f(m) = 1. For 
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Spherocylinders 




Figure 3: Solutions &(0) = 4wf(d) of the Onsager integral equation for spherocylinders for A* = 
11.01, 10.01, 9.01, 8.01. 



all three models under consideration, this amounts to finding the solution of the Onsager 
integral equation 

ln(£ f(w)) = -X J dm' f{w') u(w, w'), (28) 

where the constant £ is determined from the normalisation condition. The equation always 
has an isotropic solution f{w) = constant, and may have an anisotropic solution in addition 
at certain values of A. 

We have solved the Onsager integral equation for values of A in the range of A* = 2DL 2 X G 
[8 : 12] for spherocylinders, and A* = 2ab 2 X £ [8 : 25] for ellipsoids, in steps of 0.01. 
To this end, we have applied an iterative numerical method, which is simple and reliably 



convergent |48j]. The integrals were calculated by the Gauss-Legendre quadrature, using 50 
points per integral and checking that with 100 points we get the same results. The initial 
guess for the highest value of A*, A* = 12 or A* = 25, was 



HO rh\ = 1 exp(P 2 (cos6Q) 



where P2(u) is the second Legendre polynomial. At lower A*, the iteration was started with 
the anisotropic solution for the previous, next higher, value of A*. The convergence criterion 
was 



Up 

£ 

Li=l 



[fm+l(0i) — fm(0i)Y 



< 10" 



(29) 



where n p refers to the number of points used to perform integrals and m to the mth iteration. 
The dependence on (f> has been omitted because the distribution function does not actually 
depend on it. The resulting distribution functions for A* = 11.01, 10.01, 9.01, 8.01 for 
spherocylinders are shown in Figure 3. Our orientational distribution functions compare 
well with those in reference [49|. Due to the particular shape of the covolume between a 
pair of spherocylinders, the solutions depend on A*, but not on x. We found anisotropic 
solutions for values A* > 



.88, in agreement with Lasher []50| 
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Figure 4: Solutions = 4irf(9) of the Onsager integral equation for ellipsoids for k = 4.0, 5.0, 
10.0, 20.0 at fixed value of A* = 15.01. 



In the case of ellipsoids, the solution for fixed A* depend on the shape anisotropy parameter 
k. In Figure 4, solutions at different value of k are shown for A* = 15.01. The lowest value 
of A* which still yields an anisotropic solution is a decreasing function of k. 

Next we must identify the stable phases and the coexistence line. At coexistence, both the 
pressure and the chemical potential are equal: 

Pn = Pi] Pn = Pi- (30) 

We have solved these two equations with Newton's method, i.e., we found the density pf 
at which an isotropic phase has the pressure Pn and the density p^ at which an isotropic 
phase has the chemical potential equal to pn- At the coexistence line, both densities must 
be equal. Hence we calculated the nematic pn at which \pf — p^\ was minimal. Usually the 
value of the minimum was in the range \pf — p^\ 6 [10~ 9 ; 10~ 5 ]. 

The results are shown as a function of the anisotropy in Figure 5 for spherocylinders and in 



Figure 6 for ellipsoids. They are compared with simulation data of Bolhuis and Frenkel |38[, 



Frenkel and Mulder g9|, and Camp et al. [§2| 



The figures demonstrate that our functional performs much better than the Onsager func- 
tional, and only slightly worse than the Parsons-Lee functional. As has already been demon- 
strated elsewhere for spherocylinders and ellipsoids 0, the predictions of the latter are 
almost exact. The slightly inferior performance of our functional is probably related to 
the fact that the Parsons-Lee approach is based on the Carnahan-Starling equation of state, 
whereas the Rosenfeld functional yields the Percus-Yevick equation of state for hard spheres, 
which is slightly inferior. On principle, one can modify the Rosenfeld functional such that it 
reproduces the Carnahan-Starling equation of state for fluids, along the lines of an approach 
suggested by Tarazona [j5l| . Unfortunately, this is done at the expense of a less accurate 
description of crystals R51|| . 

Our approach could possibly be improved by making the coefficient of $3 in Eqn. ([Til) 
dependent on the orientation distribution f{w): formally, the Rosenfeld functional has 
the form of a third order y-expansion fl5^| , Mulder and Frenkel [|^| and Tjipto-Margo 
and Evans |l^] have generalized the latter to convex anisotropic bodies and applied it to 
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Figure 5: Isotropic-Nematic transition line of hard spherocylinders as a function of elongation D/L 
and reduced density p* = p/p C p (Pep is the density of a close-packed lattice). Solid lines show 
prediction of our approach j25| ) , dotted line the Onsager result @ , and dashed line the Parsons-Lee 
prediction (^). Filled squares are simulation results from Bolhuis and Frenkel |3^| (for larger D/L), 
only one transition density is given in this reference). 



Ellipsoids 
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Figure 6: Isotropic-Nematic transition line of hard ellipsoids as a function of elongation a/b and 
reduced density p* = p/p C p (Pep is the density of a close-packed lattice). Solid lines show prediction 
of our approach , dotted line the Onsager result (||) , and dashed line the Parsons-Lee prediction 
(0). Filled squares are simulation results from Camp et al. |^2|. The points at k = 3.0 are taken 
from Ref. @ 
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hard ellipsoids, determining the third virial coefficient numerically. Their nematic-isotropic 
coexistence densities were lower than those observed in the simulations. In contrast, our 
functional tends to overestimate the coexistence densities. This suggests that introducing 
orientation dependent coefficients will shift the transition lines in the correct direction. 

However, these modifications complicate the functional. We have shown that already our 
simple version describes uniform fluids reasonably well. Compared to the Parsons-Lee the- 
ory, it has the advantage of being based on a reference density functional which describes 
the local structure of hard spheres very accurately. It interpolates on the level of the direct 
correlation function, i. e. the local bulk structure, between the Percus-Yevick solution for 
hard spheres, which is very good, and the Onsager solution for infinitely elongated parti- 
cles, which is exact. Therefore, we believe that it will be suited to describe inhomogeneous 
anisotropic fluids. 



5 Summary and Outlook 



We have introduced a new density functional for liquid crystals, which interpolates between 
the successful Rosenfeld functional for hard sphere fluids |2| and the Onsager functional 
for liquid crystals [12|. As a first test, we have calculated the nematic-isotropic phase 
diagram for hard spherocylinders and hard ellipsoids and shown that the new functional 
produces reasonable results. In the next step, it shall be applied to calculate the local liquid 
structure in homogeneous fluids. For example, direct correlation functions in nematic fluids 
of ellipsoids have been determined recently from computer simulations [54, 55, o6|]. They can 
also be calculated from our density functional. This will be a second, much more sensitive 
test of the theory. 

If the new functional is successful, it will allow to predict the structure and the properties 
of nematic liquid crystals in the vicinity of inhomogeneities. It will thus contribute to 
an improved microscopic understanding of surface phenomena such as surface anchoring 
on rough and structured surfaces, or of defect structures and defect interactions. It also 
provides a new approach to smectic ordering, and might lead to useful insight into the 
nature of the nematic-smectic transition 
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